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Abstract 

>. 

The three-body decay B°(B S ) — > r yuv can occur via penguin and box diagrams in 

0\ ■ 

the standard model (SM). These channels are useful to determine the decay constants 
Jb {Jbs) and B (B s ) meson wave function. Using the B meson wave function deter- 
mined in hadronic B{B S ) decays, we calculate and get the branching ratio of order 10~ 9 
and 1CP 8 for B° and B s decay, respectively. They agree with previous calculations. 



X 



1 Introduction 

The flavor changing neutral current process is one of the most important field for testing the 
Standard Model (SM) at loop level and for establishing new physics beyond that. The rare 
B decays provide a direct and reliable tool for extracting information about the fundamental 
parameters of the Standard Model (SM), such as, the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements Vu and Vt s , if we know the value of the decay constant fs from other 
methods. Conversely, we can determine the decay constant fs if the CKM matrix elements 
are known. 



*This work is partly supported by National Science Foundation of China. 
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Pure leptonic decays B s — > and B s — > e + e~, are difficult to be measured in 

experiments, since helicity suppression give a very small branching ratio at order O(10~ 9 ) 
and (9(lCr 14 ), respectively pQ. For B° meson case the situation gets even worse due to 
the smaller CKM matrix elements V t d- For decay B s — > t + t~, although its branching 
ratio is about 10~ 7 2], it is still hard for experiments due to the low efficiency of r lepton 
measurements. 

The B°(B S ) —>■ vv decay is forbidden due to massless neutrino. Fortunately, having 
an extra real photon emitted, the radiative leptonic decays can escape from the helicity 
suppression, so that larger branching ratio of B°(B S ) — > ^vv is expected. A preliminary 
work of this type decay was carried out with many different approaches both in SM [SHU E] 
and beyond SM jBJ. In above work, it was shown that the diagrams with photon radiation 
from light quarks give the dominant contribution to the decay amplitude, that is inversely 
proportional to the constituent light quark mass. However the "constituent quark mass" 
is poorly understood. In this work, we calculate the branching ratio using B meson wave 
function which describes the constituent quark momentum distribution. The wave function 
has been studied for many years |7 4 and used in calculating non-leptonic of B decay jH]. 
Recently, this approach is also used to calculate radiative leptonic decay of charged B meson 

m 

In the next section we analyze the relevant effective Hamiltonian for the B°(B S ) — > ^vv 
decay. In SectionJSl we give our analytical and numerical results, and then compare with 
other results. At last, we summarize this work in section 0] 

2 Effective Hamiltonian 

Let us first look at the quark level process b — > qw, with q = s or d, which is shown in Fig. ^ 
This is a flavor changing neutral current process, and both box and Z penguin contribute 
to this process. The effective Hamiltonian in SM is given [TU] : 

H = C{q lll P L b){V 1 »P L v), (1) 
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Figure 1: Leading order Feynman diagrams in SM for b — > qvv, with q = s or d. 



with P L — (1 — js)/2. The coefficient C is 

(2) 
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and x = rn^/m^y. From this expression, we can see that the coefficient C is sensitive to the 
mass of the particle in loop. If new particle exist, it should affect the Wilson Coefficient 
and change the branching ratio. That is why this kind of flavor changing neutral current 
processes is sensitive to new physics 0. 

We have already mentioned that the pure leptonic (vi>) decay is forbidden due to helicity 
conservation. However, when a photon is emitted from any charged line of b or q quark, 
this pure leptonic processes turn into radiative ones and helicity suppression does not exist 
anymore. At quark level the process B s u\ — > ^vv is described by the same diagrams as 
b — > q^vv shown in FigJ21 Incidentally, we should note the following peculiarities of this 
process: 

• when photon emitted from internal charged particles (W or top quark), the above 
mentioned process will be suppressed by a factor m^jm^ (see j3j), in comparison to 
the process b —>■ quu, one can neglect the contribution of such diagrams. 

• The Wilson coefficient C is the same for the processes b — > q'yvv and b —* qvv as 
a consequence of the extension of the Low's low energy theorem (for more detail see 





Figure 2: Feynman diagrams for b — > q^vv using effective four fermi operators. 



So, when photon emitted from initial b or light quark line, there are only two diagrams, 
contributing to the process b — > q^vv. From FigHJ the corresponding decay amplitude turns 
out to be 



A 
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(3) 



3 Analytical and Numerical results 



In order to calculate analytic formulas of the decay amplitude, we use the wave functions 



decomposed in terms of spin structure. In the summation procedures, the B meson 



is treated as a heavy-light system. Thus, the B meson light-cone matrix element can be 
decomposed as |T2"] : 



B,aP 



A iP 



(4) 



where N c = 3 is color degree of freedom, Pg is the corresponding meson's momentum, 0^ 
and 0g are Lorentz scalar distribution amplitudes. As heavy quark effective theory leads to 
0^ ~ Mb<Pb, then B meson's wave function can be expressed by 



®B,a/3{x) 



[Pb + M b ] 7 5a /30i?(x). 



(5) 



In above function, the function <pB describes the momentum distribution amplitude. Since 
b quark is much heavier than the light quark in B meson, there is a sharp peak at the small 
x region for the light quark momentum fraction, 



i>_g(x) = Nbx 2 (1 — x) 2 exp 



M 2 x 2 



(6) 



It satisfies the normalization relation: 



o 



d>B(x)dx = — jB , (7) 



with J'b the B meson decay constant. This choice of B meson's wave function is almost a 
best fit from the B meson non-leptonic two body decays jH]. 

For simplicity, we consider the B meson at rest and use the light-cone coordinate (p + , p~ ,p±) 
to describe the meson and quark's momenta, where p ± = 4g (p° ±p 3 ) and p± = (p 1 ,p 2 ). Us- 
ing these coordinate we can take the B q , vv (momentum sum) and photon's momenta as 

M B , - v M B , 9 - , M B , 9 - s 

P B = -|(1,1,0 ± ); P, p = - ? |(l,r 2 ,0 ± ); P 7 = -|(0, 1 - r 2 , ± ), (8) 

with r 2 = P 2 P /M B . The momenta of b and q quark in B meson are pb = (1 — x)Pb, Pq = xPb- 
Using above convention, the amplitude for B q — > ^vv decay is written by: 

A = [ iCitvp^pPPZ + C 2 (P 7M < - P^e*) P V B ] fa-fpiMi) (9) 

with 

JO Pb'Pj Jo Pq-P~i 

Jo Pq'i r i Jo Pb ' r i 

After squaring the amplitude and performing the phase space integration over one of the 

two Dalitz variables, and summing over three generation of neutrinos, we get the differential 

decay width versus the photon energy P 7 : 

dT 6C 2 a 
dRy ~ (12vr) 2 

with 



(C[ 2 + C'i) (M B - 2P 7 )P 7 , (12) 



cs^f^-fm^. (14) 



By integrating the variable E~, we get the decay width 



r = W (cf + c?) ' (15) 
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In this work, we use following parameters EI] : 

uj b = 0.4, f B = 0.19 GeV, r B o = 1.54 x 10~ 12 s; 
uJ bs = 0.5, f Bs = 0.24 GeV, r Ba = 1.46 x lO" 12 *; 
G F = 1.66 x 10~ 5 GeV" 2 ; sin 2 U = 0.23; a = — 

V tb = 0.999, V td = 0.0074, 7 ta = 0.041. (16) 

Using these parameters, we get the branching ratios: 

Brf^B -> 7 z/z/) = 0.7 x 10~ 9 , 

Br(5 s ^ 7 w) = 2.4xl0~ 8 . (17) 

Just as we mentioned above, many approaches have been used to analyze these processes 
such as constituent quark model [3], pole model [HUE], QCD sum rule [I], and light front 
approach [S]. Here we compare our results with them which is shown in Table. ^ 



Table 1: Comparison of results from different approaches 



Mode 


Our Results 


Quark Model 


Pole Model 


Sum Rule 


light front 


BR(B° -> juu) 


0.74 x 10~ 9 


1.7 x 10~ 9 


2.1 x 10~ 9 


4.2 x 10~ 9 


1.4 x 10~ 9 


BR(B S -> 7z/z/) 


2.4 x 10~ 8 


3.5 x lO" 8 


1.8 x 10~ 8 


7.5 x lO" 8 


2.0 x lO" 8 



In constituent quark model, the non-relativity character is considered. If we replace our B 
meson distribution amplitude in eq.(JH]) by a 5 function (f B /2y/Q) 5(x — m q /m B ), our formula 
will return back to the constituent quark model in ref.jH]. Since we have poor knowledge 
about quark mass up to now, our new calculation is surely an improvement. In paper [3], 
the authors also calculate these processes in the pole model, and the results are similar to 
the quark model case. From the table, one can also see that most of the methods get similar 
results except the QCD sum rule approach, whose result is larger than others. We hope the 
experiments in future can test these different methods. 

In Fig. and Fig. (J3J), we figure out the differential decay rate of B°(Bs) — > ^vv versus 
photon energy E 1 . We also display the photon energy spectrum from constituent quark 
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Figure 3: Differential decay rate of B° — ► 71/1/ versus the photon energy £' 7 . 

model 1 . From these figures, we find our results are smaller than the constituent quark one, 
but the shape of the spectrum is the same. If normalized decay rate is used, the two lines 
will become only one, since the function is very simple 



which can be extracted from eq. (|12l) . 

Of course, there are also uncertainties in our calculation. The most large uncertainty 
comes from the heavy meson wave function. The high order twists contribution, the high 
Fock states are also not included, because they are not clear completely now. We hope that 
the non-leptonic B meson decay can offer more information in the near future. 

4 Summary 

In this work, we calculate the branching ratios in SM for B s — > jut? to be 10~ 8 and for 
B° — > ji/u to be 10~ 9 using B meson wave function constrained from non-leptonic B decays. 
1 The results of constituent quark model are from updated parameters. 
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Figure 4: Differential decay rate of B s — > ^vv versus the photon energy E 1 . 

These decay channels are useful to determine the decay constants fs and wave function. 
After calculation, we find our leading order results are at the same order as other approaches 
but a little smaller. These rare decays are sensitive to any new physics contributions which 
can be measured by future experiment such as LHCb. 
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